Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  847 1824 (N63E25 L= 13)  1397    163        1 ( 786/1801- 926/1848)   6.3 E36-E13 N59-N65
      2 1176 1733 (N53W15 L= 53)   502     79        2 (1129/1710-1216/1750)  11.3 W10-W20 N51-N54
      3  541 1473 (N32E38 L=  0)   426    118        2 ( 503/1431- 578/1528) -37.9 E43-E33 N28-N35
      4 1182 1391 (N27W11 L= 49)   736    163        4 (1121/1331-1230/1435) -49.4 W06-W13 N22-N29
      5  840 1309 (N21E12 L= 26)   269    130        1 ( 797/1282- 893/1338) -15.9 E15-E08 N19-N23
      6 1762  855 (S08W53 L= 91)   924    131        1 (1721/ 828-1828/ 886)   5.2 W49-W60 S09-S06
      7 1410  751 (S14W25 L= 63)   858    300        0 (1358/ 695-1489/ 834) -71.6 W21-W30 S17-S08
      8  822  687 (S18E13 L= 25)   775    161        1 ( 750/ 658- 899/ 715)  15.3 E18-E07 S19-S15
      9  380  525 (S31E53 L=344)  2905    332        9 ( 292/ 436- 477/ 630) -55.3 E60-E44 S37-S22
     10 1394  569 (S26W26 L= 64)  1090    295        0 (1282/ 536-1531/ 600)   3.0 W18-W37 S28-S23
  

    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