Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1092 1729 (N39E20 L= 96)  2361    470        4 ( 952/1584-1194/1829) -17.7 E34-E09 N29-N43
      2  469 1605 (N14E55 L= 61)  3799    701       12 ( 286/1439- 617/1740)  56.9 E66-E46 N00-N24
      3  643 1603 (N17E41 L= 75)   321     59        1 ( 609/1585- 673/1625)  63.4 E42-E39 N15-N19
      4  709 1471 (N11E31 L= 85)  1736    521        9 ( 597/1353- 831/1555)  14.0 E40-E24 N04-N16
      5 1529 1234 (N21W24 L=140)   716    114        2 (1481/1195-1571/1285) -23.2 W20-W27 N16-N23
      6  277 1180 (S15E53 L= 63)   417    131        0 ( 241/1123- 304/1243)  77.5 E55-E52 S19-S10
      7 1616  436 (S19W60 L=176)   499     96        1 (1576/ 417-1658/ 459)  45.0 W57-W63 S21-S16
      8 1243  238 (S44W49 L=165)  1264    184        3 (1176/ 216-1330/ 275)  31.0 W42-W54 S46-S39
  

    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