Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  756 1859 (N72E36 L=151)   963     84        2 ( 735/1818- 783/1890) -26.6 E54-E25 N66-N77
      2  405 1486 (N38E50 L=136)   340     52        1 ( 380/1463- 426/1508)  26.6 E53-E47 N36-N39
      3 1146 1471 (N32W13 L=199)   344    118        0 (1084/1461-1207/1482) -12.5 W08-W17 N31-N33
      4  646 1465 (N36E25 L=161)   299     72        0 ( 609/1447- 682/1481) -18.4 E28-E22 N34-N36
      5 1049 1367 (N26W04 L=191)   631    176        1 (1014/1286-1080/1446) -84.3 W02-W06 N21-N31
      6  343 1367 (N30E52 L=135)   391     81        1 ( 334/1324- 357/1408) -78.7 E54-E50 N27-N32
      7 1442 1373 (N23W32 L=219)   359    156        1 (1378/1334-1501/1402) -40.6 W27-W35 N19-N25
      8 1752 1217 (N09W54 L=241)   907    154        3 (1729/1152-1782/1287) -90.0 W52-W56 N05-N13
      9  699 1171 (N16E20 L=167)   754    250        1 ( 643/1088- 755/1243) -90.0 E23-E16 N11-N20
     10 1750  857 (S13W50 L=237)   724     80        2 (1711/ 829-1781/ 885)   8.1 W46-W54 S14-S10
     11 1534  636 (S24W33 L=219)  1426    290        1 (1458/ 507-1584/ 747)  65.0 W28-W36 S32-S17
     12 1718  708 (S22W49 L=236)   455     97        0 (1674/ 676-1754/ 744)  20.6 W44-W52 S23-S19
     13 1767  668 (S26W56 L=242)  1218    206        0 (1694/ 601-1832/ 735)  26.6 W49-W63 S29-S21
     14  227  647 (S16E69 L=117)  1702    143        7 ( 182/ 599- 278/ 713) -36.9 E76-E62 S19-S11
     15  881  595 (S21E13 L=174)   483    155        0 ( 811/ 554- 948/ 634)  26.6 E18-E07 S23-S18
     16  967  442 (S33E09 L=177)   300     72        0 ( 949/ 402- 983/ 477)  51.3 E11-E07 S35-S30
     17  875  351 (S40E18 L=168)   305     80        1 ( 842/ 332- 918/ 363)   8.1 E21-E14 S41-S38
  

    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