Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  967 1842 (N43E37 L=145)   540     64        0 ( 931/1824-1006/1859)   0.0 E41-E32 N41-N44
      2  615 1664 (N20E47 L=135)   345     82        1 ( 577/1641- 645/1696)  68.2 E48-E45 N17-N22
      3  934 1559 (N21E21 L=161)   554    270        0 ( 900/1461- 987/1660) -40.2 E27-E14 N16-N27
      4 1320 1583 (N35W01 L=183)  1141    213        3 (1252/1509-1394/1640) -17.1 E05-W08 N32-N37
      5  915 1409 (N12E17 L=165)   404    114        0 ( 868/1368- 959/1457) -23.2 E20-E13 N10-N13
      6 1709 1170 (N22W39 L=221)  1361    342        1 (1670/1017-1732/1321) -64.8 W33-W43 N13-N30
      7  572 1261 (S05E33 L=149)   285     83        0 ( 551/1222- 592/1301) -78.7 E33-E31 S07-S02
      8  681 1191 (S07E24 L=158)   737    275        0 ( 581/1111- 736/1293)  67.4 E28-E22 S12-N00
      9 1837 1224 (N31W52 L=234)   903    152        0 (1784/1159-1865/1272) -24.0 W45-W55 N27-N33
     10 1199 1166 (N06W06 L=188)   646    153        0 (1142/1110-1267/1201)  63.4 W03-W08 N01-N09
     11 1421 1112 (N09W20 L=202)   484     86        0 (1397/1068-1449/1156)  80.5 W18-W21 N06-N12
     12  180 1029 (S26E64 L=119)  1138    130        1 ( 131/ 970- 229/1066)  21.8 E72-E56 S30-S23
     13  251  825 (S38E54 L=128)   773    142        2 ( 201/ 793- 315/ 865)  17.4 E61-E45 S40-S34
     14  865  842 (S21E04 L=178)   212     67        0 ( 826/ 826- 895/ 857)  21.8 E06-E01 S22-S19
     15 1150  585 (S27W22 L=204)   381     70        1 (1122/ 556-1181/ 612)   0.0 W19-W24 S28-S25
     16 1501  294 (S30W68 L=250)  1478    172        2 (1425/ 268-1567/ 324)  13.4 W59-W81 S33-S27
  

    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