Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1510 1742 (N31W59 L=202)  2305    448        5 (1407/1688-1613/1777) -10.8 W47-W69 N26-N32
      2 1040 1054 (N06W02 L=145)   319    108        0 (1005/ 999-1095/1099)  26.6 E01-W05 N03-N07
      3  339 1012 (N20E44 L=100)   578    158        1 ( 314/ 934- 374/1087)  53.1 E46-E39 N14-N23
      4 1196 1029 (N00W09 L=153)   288     98        0 (1148/ 984-1241/1058) -45.0 W07-W11 S03-N02
      5 1840  898 (S26W52 L=195)   701    152        0 (1806/ 824-1872/ 967)  41.2 W47-W56 S29-S22
      6  971  838 (S05E08 L=136)   181     58        0 ( 947/ 807- 989/ 861) -76.0 E08-E06 S06-S02
      7 1595  721 (S28W28 L=171)   660    207        0 (1534/ 623-1651/ 831)  39.3 W21-W33 S33-S23
      8 1168  502 (S30E06 L=138)   297     57        0 (1135/ 482-1198/ 527)  11.3 E08-E03 S30-S28
      9  489  364 (S21E57 L= 87)  2383    536        3 ( 445/ 238- 521/ 495)  49.9 E66-E48 S31-S12
     10  569  437 (S18E46 L= 98)   296     46        0 ( 550/ 414- 586/ 459) -90.0 E47-E43 S19-S16
     11  650  348 (S26E45 L= 99)   635    161        0 ( 634/ 292- 667/ 409)  48.8 E48-E41 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