Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1071 1789 (N42E26 L=200)   568    180        0 (1026/1703-1118/1851) -27.9 E34-E17 N37-N46
      2 1013 1443 (N17E13 L=213)   792    220        0 ( 937/1355-1084/1499) -18.4 E18-E06 N13-N18
      3  234 1467 (N00E74 L=152)   709     75        0 ( 216/1431- 260/1496) -56.3 E77-E69 S01-N02
      4 1864 1049 (N22W57 L=283)  1160    181        0 (1806/ 966-1896/1100) -18.4 W49-W63 N17-N24
      5  619  995 (S19E23 L=203)   900    369        1 ( 466/ 890- 763/1044)  45.0 E31-E14 S28-S12
      6  279  836 (S37E50 L=176)   424     80        0 ( 263/ 793- 303/ 876)  56.3 E51-E46 S39-S33
      7  715  823 (S26E13 L=213)   271     87        0 ( 679/ 796- 753/ 844)  51.3 E14-E10 S28-S23
      8  824  687 (S31E02 L=224)  1320    327        3 ( 671/ 659- 956/ 717)  27.9 E12-W05 S35-S26
      9 1082  603 (S28W17 L=243)   510    162        1 (1010/ 563-1152/ 651)   0.0 W11-W22 S29-S26
     10 1454  531 (S20W43 L=269)   550    129        0 (1399/ 487-1515/ 574)  76.0 W40-W45 S23-S15
     11  346  479 (S62E53 L=173)   514     43        2 ( 324/ 463- 370/ 505) -33.7 E61-E46 S63-S59
     12 1534  324 (S27W67 L=293)   585     94        0 (1492/ 292-1573/ 358)  80.5 W65-W68 S30-S24
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y