Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  488 1659 (N19E57 L=274)  4223    704       11 ( 431/1495- 558/1789) -49.6 E65-E46 N09-N29
      2 1592 1598 (N50W26 L=356)   764    120        1 (1555/1566-1627/1631)  41.6 W21-W30 N45-N53
      3  939 1493 (N20E18 L=312)   279    109        0 ( 883/1471- 967/1532)  45.0 E21-E16 N17-N22
      4 1340 1244 (N16W12 L=343)   584    143        2 (1299/1170-1390/1317) -36.9 W08-W16 N13-N19
      5 1849 1164 (N28W53 L= 24)   469     68        2 (1830/1134-1867/1196) -36.9 W50-W55 N26-N29
      6  230 1135 (S17E57 L=274)  1152    457        5 ( 110/1105- 326/1160)  15.3 E71-E48 S20-S14
      7 1337  871 (S05W22 L=352)   221     71        1 (1301/ 849-1362/ 905)  76.0 W20-W23 S06-S02
      8 1404  631 (S16W34 L=  5)   233     81        0 (1360/ 616-1442/ 648)  45.0 W31-W36 S17-S13
      9 1551  358 (S26W60 L= 31)  1494    266        1 (1460/ 337-1633/ 386)  23.2 W52-W67 S29-S21
  

    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