Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  604 1760 (N58E29 L=160)   919    176        0 ( 532/1722- 678/1800)  12.5 E38-E20 N55-N60
      2 1159 1702 (N41W23 L=213)   704    119        2 (1105/1681-1215/1731) -12.5 W18-W27 N39-N44
      3 1319 1623 (N33W32 L=221)   611     83        0 (1262/1606-1366/1644)  -7.1 W27-W35 N31-N34
      4 1665 1536 (N22W56 L=246)  1153    151        3 (1600/1496-1739/1571) -35.0 W50-W62 N18-N25
      5 1418  897 (S13W22 L=212)   249    117        1 (1365/ 877-1479/ 918)  -7.1 W18-W26 S13-S11
      6 1021  726 (S18E04 L=185)   896    362        0 ( 949/ 589-1065/ 883)  87.0 E07-E01 S26-S07
      7 1653  794 (S22W39 L=229)   528    175        0 (1615/ 720-1678/ 874)  60.9 W36-W41 S26-S17
      8 1442  602 (S32W23 L=212)   437     77        0 (1403/ 583-1478/ 631)   9.5 W19-W25 S33-S30
      9  418  552 (S20E51 L=139)   437     99        0 ( 406/ 495- 439/ 593)  45.0 E53-E47 S23-S17
     10 1482  383 (S49W29 L=219)   985    140        0 (1414/ 360-1541/ 409) -14.9 W21-W37 S52-S46
  

    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