Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1020 1905 (N52E44 L= 90)  1223    156        0 ( 966/1857-1051/1943) -31.4 E56-E36 N46-N57
      2 1082 1835 (N47E30 L=105)   577     60        1 (1060/1799-1099/1867) -29.1 E34-E25 N44-N50
      3  697 1771 (N30E50 L= 84)   591    117        1 ( 654/1715- 746/1819) -19.4 E59-E42 N26-N32
      4 1432  916 (S00W26 L=160)  1188    359        0 (1344/ 789-1492/1066) -27.9 W16-W33 S05-N05
      5  538  855 (S29E26 L=109)   265     61        1 ( 500/ 833- 571/ 877)  31.0 E28-E23 S30-S26
      6  382  432 (S65E49 L= 85)   536     53        1 ( 355/ 411- 400/ 461)  15.5 E59-E41 S66-S61
      7 1502  365 (S26W59 L=194)   662    146        0 (1438/ 338-1571/ 388)  45.0 W55-W63 S29-S22
      8 1281  343 (S36W45 L=179)   822     80        1 (1228/ 312-1325/ 376)  51.3 W42-W46 S38-S32
      9  881  207 (S62W34 L=168)   582     56        0 ( 848/ 183- 909/ 236) -80.5 W30-W39 S64-S58
  

    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