Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1562 1718 (N57W54 L=162)   801     96        1 (1532/1691-1591/1741)  15.6 W44-W69 N53-N60
      2 1456 1454 (N33W27 L=134)   850    292        0 (1329/1379-1543/1527) -25.2 W16-W33 N28-N37
      3 1005 1389 (N24E06 L=102)   946    401        1 ( 915/1277-1127/1503) -39.1 E13-W03 N17-N30
      4  913 1316 (N18E10 L= 97)   500    150        1 ( 895/1236- 935/1374) -71.6 E12-E08 N12-N21
      5 1127 1312 (N20W03 L=111)   679    164        1 (1060/1249-1195/1343) -26.6 E01-W08 N16-N21
      6  832  761 (S17E09 L= 98)   578    117        1 ( 784/ 723- 896/ 789)  29.7 E12-E05 S19-S14
      7 1386  609 (S21W29 L=137)   488    123        1 (1357/ 572-1420/ 650) -18.4 W26-W32 S23-S18
      8  258  598 (S36E65 L= 43)   686     95        0 ( 222/ 566- 305/ 618)  13.5 E80-E55 S39-S33
      9 1439  389 (S36W43 L=150)   653    114        1 (1391/ 355-1488/ 427)  56.3 W40-W45 S38-S32
  

    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