Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1155 1487 (N26W14 L= 49)  1023    238        0 (1030/1441-1256/1514) -23.2 W05-W19 N22-N28
      2  447 1366 (N26E37 L=359)  1467    328        4 ( 319/1281- 568/1456) -32.0 E50-E26 N18-N33
      3 1003 1155 (N07E00 L= 35)   258     83        0 ( 969/1117-1036/1198) -68.2 E01-W01 N04-N09
      4 1241 1109 (N02W14 L= 49)   280    108        1 (1199/1069-1271/1158)  33.7 W10-W16 N00-N04
      5  764 1047 (N03E15 L= 20)   303    140        0 ( 742/ 975- 785/1109)  74.1 E16-E13 S01-N06
      6 1847 1017 (S09W59 L= 94)   827    180        1 (1810/ 937-1870/1099)  52.1 W54-W62 S13-S04
      7 1907  913 (S16W68 L=103)  1315    109        3 (1876/ 879-1942/ 954) -45.0 W62-W74 S17-S12
      8  615  517 (S28E35 L=  1)  1319    261        1 ( 548/ 411- 681/ 633)  37.9 E45-E27 S34-S20
  

    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