Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  479 1646 (N18E56 L=275)  4311    718        5 ( 424/1487- 547/1783) -49.6 E66-E43 N08-N28
      2 1580 1598 (N49W24 L=356)   776    141        2 (1533/1562-1618/1628)  38.7 W19-W29 N45-N53
      3  919 1493 (N19E19 L=312)   215     93        0 ( 864/1478- 951/1517)  36.9 E22-E17 N17-N21
      4 1336 1235 (N15W12 L=344)   789    269        1 (1278/1093-1428/1349) -35.5 W06-W21 N10-N20
      5 1843 1172 (N28W52 L= 24)   527     75        1 (1820/1144-1863/1203) -33.7 W49-W54 N26-N29
      6  225 1142 (S17E58 L=274)  1571    473        2 ( 105/1103- 322/1166)  16.9 E71-E48 S20-S13
      7 1396  641 (S15W33 L=  5)   231    103        0 (1349/ 625-1443/ 664)  45.0 W30-W35 S17-S12
      8 1687  585 (S09W57 L= 29)   347    105        0 (1653/ 560-1706/ 612)  76.0 W55-W59 S11-S07
      9 1548  364 (S26W59 L= 31)  1533    276        1 (1444/ 347-1636/ 403)  34.7 W51-W66 S29-S20
  

    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