Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1097 1787 (N43E21 L=172)   966    225        1 (1054/1719-1146/1859) -33.7 E28-E13 N38-N48
      2  536 1608 (N15E48 L=145)   545    178        2 ( 520/1542- 550/1664) -45.0 E52-E43 N11-N19
      3 1226 1115 (N03W09 L=202)   278     84        1 (1182/1096-1274/1129)  31.0 W06-W11 N01-N04
      4 1388 1046 (N04W20 L=213)   322     56        1 (1364/1011-1413/1074)  76.0 W17-W21 N01-N05
      5  229 1027 (S24E54 L=138)   635    205        1 ( 171/ 977- 289/1070)  35.0 E59-E49 S27-S20
      6 1414  866 (S05W26 L=219)   910    165        1 (1384/ 781-1433/ 953) -53.1 W23-W29 S09-N00
      7 1181  807 (S15W15 L=207)  1299    282        2 (1036/ 750-1313/ 875)  45.0 W07-W20 S21-S07
      8 1844  821 (N06W61 L=254)   393    140        1 (1807/ 804-1881/ 837)  11.3 W56-W66 N04-N07
      9  858  532 (S40W04 L=197)  1701    343        3 ( 692/ 503-1035/ 566)  15.6 E09-W16 S44-S35
     10 1450  468 (S25W45 L=237)  1274    249        0 (1327/ 443-1549/ 500)  32.7 W36-W51 S29-S20
     11 1572  413 (S24W57 L=250)   405     24        0 (1552/ 390-1593/ 437) -71.6 W55-W58 S25-S22
  

    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