Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  966 1747 (N49W10 L=164)   314     14        0 ( 951/1723- 985/1768)  63.4 W08-W11 N46-N50
      2  844 1562 (N36E04 L=150)   257     78        0 ( 802/1550- 885/1572) -31.0 E06-E01 N34-N37
      3  733 1540 (N36E12 L=142)   310     80        0 ( 701/1502- 767/1571)  33.7 E15-E09 N34-N38
      4 1297  727 (S21W13 L=167)  1126    388        0 (1122/ 662-1431/ 825)   5.4 W01-W22 S24-S17
      5  403  643 (S14E48 L=106)   547     66        1 ( 370/ 607- 426/ 682)  38.7 E51-E45 S15-S11
      6  840  638 (S20E17 L=137)   208     35        0 ( 822/ 616- 865/ 658)  33.7 E18-E15 S21-S19
      7  215  554 (S17E78 L= 77)  1551     78        0 ( 193/ 509- 234/ 590)  32.0 E84-E72 S19-S14
      8  536  385 (S32E50 L=104)  1689    251        0 ( 459/ 335- 612/ 448) -33.7 E56-E41 S35-S26
  

    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