Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  469 1708 (N70E47 L=108)   550     22        2 ( 450/1688- 490/1729) -45.0 E57-E39 N68-N73
      2 1443 1579 (N26W46 L=201)  1608    334        2 (1310/1541-1595/1609) -22.2 W35-W58 N20-N30
      3 1579 1471 (N16W50 L=205)   482    124        1 (1533/1434-1616/1494) -56.3 W46-W52 N12-N18
      4 1590 1388 (N11W46 L=201)   362    104        0 (1548/1358-1634/1416) -51.3 W43-W48 N08-N13
      5  480 1261 (N34E31 L=124)   327     80        0 ( 438/1235- 520/1286) -35.5 E34-E27 N31-N36
      6  425 1139 (N27E37 L=118)  1006    309        1 ( 384/1044- 510/1274) -69.4 E40-E31 N20-N36
      7  605 1118 (N22E23 L=132)   227     67        1 ( 577/1094- 637/1144) -53.1 E24-E21 N20-N24
      8 1774 1055 (S13W50 L=205)   405     45        0 (1757/1023-1792/1087) -76.0 W48-W50 S14-S10
      9  338  854 (N12E49 L=106)   362    125        1 ( 297/ 825- 364/ 903)  23.2 E53-E46 N10-N14
     10  224  832 (N12E62 L= 93)   744    101        0 ( 186/ 813- 276/ 861) -16.7 E66-E56 N10-N14
     11  384  760 (N06E48 L=107)   543    142        1 ( 355/ 690- 410/ 818)  82.9 E50-E46 N01-N09
     12  644  704 (S02E31 L=125)   355     77        1 ( 615/ 671- 669/ 746) -76.0 E31-E28 S04-N00
     13  514  535 (S10E46 L=109)   309     67        1 ( 480/ 516- 544/ 560)   0.0 E49-E43 S10-S08
     14  384  460 (S13E63 L= 92)   639     66        2 ( 365/ 419- 405/ 492)  32.0 E67-E58 S15-S10
     15 1142  408 (S34E10 L=145)   317     59        1 (1113/ 382-1163/ 430) -76.0 E11-E08 S35-S31
     16  617  323 (S26E51 L=104)   410     50        2 ( 595/ 299- 645/ 346) -56.3 E54-E48 S27-S24
  

    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