Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  593 1714 (N31E52 L=345)  1257    121        2 ( 547/1655- 632/1761) -17.4 E60-E44 N27-N34
      2  400 1608 (N20E62 L=334)  1749    186        4 ( 325/1566- 491/1655)  33.7 E70-E55 N16-N24
      3  172 1209 (S08E68 L=329)  2041    162        6 ( 121/1153- 229/1268)  36.0 E76-E61 S11-S03
      4 1378 1180 (N14W18 L= 54)   450    128        2 (1350/1137-1401/1241) -71.6 W15-W19 N11-N17
      5  344 1183 (S06E48 L=349)   301     69        0 ( 314/1159- 385/1202)  -8.1 E50-E43 S06-S05
      6 1052 1071 (N01W01 L= 37)  1117    253        2 ( 979/ 962-1128/1166)  77.0 E02-W03 S05-N08
      7  220 1013 (S18E57 L=339)   821    186        1 ( 161/ 940- 289/1075)  35.5 E65-E50 S23-S13
      8 1310  861 (S05W20 L= 57)   680    151        1 (1251/ 818-1371/ 919)  58.0 W17-W22 S08-N00
      9 1591  880 (N02W38 L= 75)   318     52        1 (1574/ 850-1614/ 916) -71.6 W37-W40 N00-N03
     10 1785  754 (S01W59 L= 95)  1313    119        1 (1735/ 709-1822/ 809)  39.8 W54-W62 S04-N02
     11  213  777 (S33E60 L=337)   962     96        1 ( 170/ 755- 260/ 802)  -4.4 E66-E53 S34-S31
     12  615  572 (S38E18 L= 18)   246     54        0 ( 582/ 552- 645/ 594)  31.0 E21-E16 S40-S36
  

    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