Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1837 (N55E34 L=207)  1267    157        2 ( 854/1805-1004/1866)   0.0 E44-E22 N52-N57
      2 1277 1717 (N52W07 L=249)   349     29        0 (1254/1697-1303/1733)  21.8 W05-W10 N49-N53
      3  550 1667 (N33E52 L=189)  1006    159        1 ( 489/1620- 606/1717) -12.8 E65-E43 N30-N35
      4 1444 1600 (N45W24 L=265)   643     64        1 (1411/1569-1478/1626)  -8.1 W20-W27 N42-N46
      5  919 1520 (N29E16 L=225)   356    114        0 ( 868/1502- 975/1536)  23.2 E19-E12 N27-N30
      6  671 1514 (N24E33 L=208)   256     13        0 ( 651/1498- 696/1529)   0.0 E34-E31 N23-N25
      7 1729 1214 (N22W48 L=289)   953    184        0 (1659/1145-1785/1278) -26.6 W41-W53 N18-N25
      8  240  788 (S27E58 L=183)   572    107        1 ( 215/ 737- 257/ 846)  48.8 E62-E54 S30-S22
  

    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