Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  667 1517 (N33E29 L= 13)   668    162        0 ( 612/1471- 720/1565) -26.6 E34-E24 N29-N35
      2 1191 1328 (N22W10 L= 52)  1468    360        2 (1023/1227-1358/1436) -29.7 E01-W21 N16-N28
      3  300 1353 (N19E56 L=346)   872    137        0 ( 279/1296- 322/1426) -69.4 E58-E52 N16-N24
      4  704 1397 (N25E24 L= 18)   331     50        1 ( 681/1378- 737/1419)  33.7 E25-E21 N23-N26
      5  859 1301 (N19E12 L= 30)   730    214        0 ( 813/1218- 931/1400) -57.5 E15-E06 N14-N25
      6  718 1270 (N17E21 L= 21)   252     73        1 ( 688/1239- 743/1301) -45.0 E23-E18 N14-N18
      7  479 1136 (N07E36 L=  5)   859    194        1 ( 450/1047- 532/1221) -63.4 E39-E31 N02-N12
      8  569 1021 (N01E29 L= 13)   237     85        0 ( 532/1009- 609/1040)  21.8 E31-E26 N00-N02
      9  394  979 (S03E42 L=360)   390     94        2 ( 361/ 958- 436/1005)  -9.5 E44-E38 S03-S01
     10  992  503 (S31E00 L= 42)   267     70        0 ( 965/ 471-1013/ 532) -63.4 E02-W01 S33-S29
     11 1183  264 (S51W20 L= 62)   603    139        0 (1132/ 228-1222/ 297) -32.5 W13-W25 S54-S47
  

    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