Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1570 (N30E59 L=180)   970    111        1 ( 349/1533- 449/1604)  -3.4 E70-E53 N27-N31
      2  442 1497 (N26E50 L=190)   343     76        2 ( 428/1458- 456/1535) -39.8 E53-E46 N23-N28
      3 1663 1463 (N35W49 L=288)   720     72        2 (1628/1426-1689/1500) -90.0 W44-W51 N32-N37
      4  438 1229 (N09E41 L=198)   378    134        1 ( 380/1216- 513/1244)   5.2 E46-E35 N07-N09
      5 1628  808 (S07W43 L=282)   972    194        1 (1596/ 730-1648/ 898) -68.2 W39-W45 S11-S01
      6  543  831 (S15E30 L=210)   670    220        1 ( 451/ 781- 637/ 866)  -8.7 E36-E23 S16-S12
      7  840  676 (S22E09 L=231)   270    116        0 ( 821/ 626- 865/ 717) -78.7 E10-E07 S24-S18
      8 1141  591 (S25W12 L=252)   304    123        0 (1086/ 574-1191/ 622)  -7.1 W07-W15 S25-S23
      9 1367  532 (S27W29 L=269)   400     92        0 (1327/ 514-1403/ 570)  51.3 W26-W31 S28-S23
     10  382  469 (S42E55 L=185)   638     70        0 ( 367/ 430- 404/ 500)  26.6 E60-E48 S45-S39
  

    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