Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1051 1807 (N56E16 L=168)   557     45        3 (1024/1783-1074/1833)  90.0 E18-E12 N53-N59
      2  968 1800 (N53E24 L=161)   532     50        2 ( 942/1770- 991/1824) -63.4 E26-E20 N50-N56
      3 1185 1787 (N57E01 L=183)   791    121        2 (1134/1761-1242/1809)   0.0 E07-W05 N54-N59
      4  579 1579 (N29E43 L=141)   581    124        2 ( 557/1509- 609/1631) -58.0 E45-E39 N24-N32
      5  321 1482 (N19E63 L=121)   843    110        4 ( 289/1434- 350/1522) -56.3 E67-E59 N15-N21
      6 1277 1451 (N31W11 L=195)   761    256        0 (1170/1385-1372/1510) -17.4 W02-W18 N27-N33
      7 1621 1346 (N29W39 L=223)   650    166        1 (1551/1306-1692/1378)  23.2 W32-W46 N25-N31
      8  478 1029 (S07E35 L=149)   434    102        1 ( 444/ 980- 501/1069) -51.3 E38-E33 S09-S03
      9  262  885 (S18E54 L=130)   623    140        2 ( 239/ 828- 295/ 934)  54.5 E57-E50 S22-S15
     10  316  834 (S21E48 L=136)   308     65        1 ( 299/ 802- 333/ 869)  68.2 E50-E46 S23-S18
  

    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