Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  506 1508 (N10E45 L=280)  3016    621        3 ( 450/1270- 592/1729) -62.7 E51-E35 S03-N26
      2 1756 1470 (N42W43 L=  8)   507    110        1 (1728/1442-1784/1490)  36.9 W39-W47 N38-N44
      3 1638 1399 (N31W31 L=356)   277     49        0 (1605/1380-1659/1413)  14.0 W28-W32 N30-N32
      4 1288 1321 (N16W08 L=333)   377    126        0 (1259/1250-1314/1383) -50.2 W05-W11 N12-N19
      5  150 1214 (S12E69 L=256)  1170    190        2 ( 124/1180- 182/1262)  45.0 E74-E64 S14-S08
      6  230 1172 (S14E58 L=267)  1121    249        0 ( 160/1149- 303/1193)   7.6 E65-E50 S15-S11
      7 1710  975 (N08W44 L=  8)   356     82        0 (1691/ 936-1726/1018) -45.0 W41-W45 N05-N09
      8  723  940 (S19E16 L=309)   300    124        0 ( 666/ 925- 782/ 972)  35.5 E19-E12 S20-S15
      9 1403  699 (S16W31 L=356)   269     60        0 (1369/ 670-1437/ 725)  59.0 W29-W32 S18-S13
     10 1698  605 (S12W56 L= 21)   403     95        0 (1671/ 571-1712/ 633) -36.9 W53-W58 S13-S10
     11 1563  442 (S25W54 L= 19)  1815    238        3 (1502/ 407-1628/ 489)  40.6 W48-W60 S27-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