Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1359 1645 (N34W33 L=180)   315     71        1 (1333/1627-1401/1665) -21.8 W30-W36 N33-N36
      2 1635 1634 (N30W55 L=203)   567    175        0 (1610/1593-1672/1658)  24.0 W51-W60 N27-N32
      3  305 1465 (N35E50 L= 97)  1929    448        3 ( 210/1435- 396/1503) -16.9 E63-E40 N32-N39
      4 1810 1420 (N15W64 L=211)   600    108        0 (1790/1385-1839/1443)  23.2 W61-W68 N12-N16
      5  550 1258 (N18E28 L=119)   246     84        0 ( 508/1245- 589/1273) -18.4 E30-E24 N17-N19
      6  184 1170 (N17E60 L= 88)   786    232        1 ( 147/1099- 230/1224) -45.0 E64-E54 N12-N21
      7  244 1149 (N15E53 L= 94)   372    107        2 ( 222/1123- 266/1201)  63.4 E55-E50 N14-N18
      8  142 1104 (N14E65 L= 82)  2002    394        6 ( 105/1010- 181/1186) -57.5 E70-E60 N08-N19
      9  643 1123 (N09E22 L=125)   393     90        0 ( 595/1082- 677/1165) -50.2 E25-E20 N06-N12
     10  834 1021 (N01E11 L=136)   366    186        2 ( 794/ 960- 866/1075)  33.7 E14-E08 S02-N04
     11 1853 1016 (S10W59 L=206)   723    132        4 (1828/ 963-1884/1073) -90.0 W56-W62 S12-S06
     12 1756  974 (S11W49 L=196)   363     79        2 (1723/ 948-1779/1004) -11.3 W45-W50 S12-S09
     13 1385  724 (S23W20 L=167)   210     50        1 (1361/ 706-1408/ 744) -33.7 W18-W21 S24-S21
     14 1248  678 (S24W10 L=158)   385    131        1 (1181/ 660-1312/ 702) -14.0 W06-W14 S25-S22
     15  465  604 (S20E43 L=104)   275     61        1 ( 442/ 585- 488/ 632)  33.7 E45-E41 S20-S18
     16  426  449 (S29E53 L= 94)  1100    185        5 ( 405/ 380- 456/ 514)  66.0 E57-E49 S33-S24
     17 1065  405 (S41E06 L=142)   400     81        1 (1024/ 387-1104/ 437)   0.0 E08-E02 S42-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