Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1230 1738 (N54W10 L= 86)   497    107        1 (1186/1710-1279/1758) -16.7 W04-W14 N51-N55
      2  500 1489 (N25E44 L= 32)  1451    332        2 ( 441/1369- 538/1628) -50.5 E52-E37 N17-N34
      3 1337 1436 (N31W18 L= 94)  2189    575        1 (1139/1243-1480/1561) -31.6 W02-W28 N19-N38
      4 1182 1225 (N15W08 L= 84)   302    157        0 (1154/1172-1228/1292) -54.5 W05-W11 N12-N19
      5 1602 1218 (N19W38 L=114)   280    126        0 (1577/1178-1628/1253) -76.0 W35-W39 N16-N20
      6 1299  670 (S18W22 L= 98)   570    167        0 (1227/ 621-1371/ 727)  41.2 W17-W25 S21-S13
      7 1715  653 (S15W56 L=132)   573     41        1 (1688/ 622-1740/ 677) -53.1 W52-W59 S17-S13
      8 1651  580 (S21W52 L=129)   452    103        0 (1627/ 547-1681/ 624)  59.0 W50-W54 S22-S17
      9 1714  594 (S19W58 L=135)   561     98        0 (1679/ 567-1760/ 615)  20.6 W54-W63 S21-S17
     10  421  485 (S41E47 L= 29)   609     64        0 ( 397/ 464- 457/ 516) -26.6 E50-E42 S43-S38
  

    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