Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  763 1513 (N42E08 L=351)  1428    286        5 ( 678/1457- 892/1605) -43.0 E14-W01 N36-N51
      2 1031 1408 (N29W09 L=  7)   202     47        1 (1006/1390-1058/1428) -63.4 W07-W10 N27-N31
      3 1196 1348 (N22W18 L= 16)   270    103        0 (1179/1303-1216/1401)  71.6 W16-W20 N19-N25
      4  240 1294 (N36E60 L=298)  1510    153        3 ( 173/1242- 307/1334)  -7.1 E74-E50 N32-N39
      5  436 1153 (N25E37 L=321)   917    152        3 ( 366/1094- 482/1204) -35.0 E43-E33 N20-N28
      6 1770  950 (S15W49 L= 47)   488    117        1 (1745/ 896-1796/1006)  50.2 W45-W51 S17-S11
      7  590  964 (N11E28 L=330)   219     66        0 ( 552/ 948- 626/ 977)   0.0 E30-E25 N10-N11
      8  191  916 (N13E64 L=294)   453     65        2 ( 172/ 879- 212/ 944) -63.4 E66-E61 N10-N14
      9  261  776 (N04E59 L=299)   452     75        3 ( 234/ 743- 287/ 802)  31.0 E62-E55 N02-N05
     10  827  690 (S09E18 L=340)   708    240        2 ( 751/ 623- 897/ 768)  26.6 E24-E12 S11-S05
     11  510  661 (S06E40 L=318)   260     77        0 ( 480/ 629- 536/ 681) -56.3 E42-E38 S08-S04
     12  531  595 (S10E41 L=317)   284     21        0 ( 508/ 576- 553/ 615) -33.7 E42-E39 S11-S08
     13 1435  560 (S31W17 L= 15)   313     97        1 (1380/ 548-1477/ 575) -15.9 W13-W20 S32-S29
     14 1712  540 (S41W45 L= 43)   558     66        2 (1691/ 506-1732/ 573) -80.5 W41-W48 S44-S38
     15 1291  434 (S36W05 L=  3)   519     92        0 (1256/ 407-1329/ 463)  26.6 W02-W08 S38-S34
     16  950  328 (S36E23 L=336)   463     88        1 ( 926/ 289- 971/ 372)  40.6 E26-E19 S39-S32
  

    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