Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1001 1710 (N45W18 L=119)   267     70        0 ( 964/1693-1035/1728)   9.5 W14-W20 N44-N46
      2 1712 1380 (N08W54 L=155)   413     74        0 (1677/1349-1745/1409) -45.0 W51-W56 N05-N09
      3  297 1083 (N20E47 L= 55)  1789    365        2 ( 200/ 988- 413/1181) -43.2 E55-E38 N11-N27
      4  135  941 (N14E67 L= 35)  1699    323        7 (  91/ 860- 176/1002) -90.0 E72-E60 N09-N18
      5 1593  700 (S29W30 L=132)   944    202        3 (1543/ 606-1634/ 812)  56.3 W25-W34 S34-S22
      6 1155  532 (S30E03 L= 99)  1008    289        1 (1117/ 415-1198/ 676)  60.6 E08-W02 S37-S21
      7  694  650 (S13E28 L= 74)   239     66        0 ( 667/ 625- 727/ 671)  14.0 E29-E25 S13-S11
      8  407  455 (S19E57 L= 45)   527     73        1 ( 389/ 427- 427/ 486)  63.4 E59-E54 S20-S16
      9  607  451 (S23E41 L= 61)   304     60        0 ( 577/ 436- 640/ 468)   0.0 E43-E37 S23-S22
     10  592  338 (S30E48 L= 54)  1200    236        0 ( 565/ 269- 625/ 415)  36.9 E54-E41 S34-S25
  

    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