Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  382 1428 (N31E42 L=240)   322     50        0 ( 355/1406- 408/1455)  33.7 E44-E39 N30-N33
      2 1449 1371 (N15W31 L=313)   251     80        0 (1437/1326-1460/1416)  71.6 W29-W32 N12-N18
      3  670 1340 (N22E19 L=263)   386    147        0 ( 608/1300- 743/1359)  12.5 E23-E14 N20-N23
      4  293 1201 (N18E48 L=233)   415     99        0 ( 257/1186- 332/1216) -15.9 E51-E44 N16-N19
      5  155 1055 (N11E64 L=218)  1411    235        6 ( 121/ 997- 197/1112)  74.1 E68-E58 N07-N14
      6 1813  999 (S10W54 L=336)   841    215        0 (1777/ 894-1845/1098)  59.7 W50-W58 S16-S04
      7 1872  939 (S14W60 L=342)   489     92        0 (1849/ 901-1898/ 984)  39.8 W57-W64 S16-S11
      8 1618  796 (S21W37 L=318)   429    139        1 (1570/ 747-1652/ 870)  49.4 W33-W39 S23-S16
      9 1350  751 (S21W18 L=299)   907    194        3 (1280/ 676-1397/ 823)  41.6 W12-W21 S24-S16
     10  549  742 (S12E33 L=248)   782    147        2 ( 496/ 705- 601/ 775) -23.2 E37-E29 S14-S10
     11  220  697 (S11E63 L=218)   576    116        2 ( 205/ 658- 241/ 747) -90.0 E65-E60 S13-S08
     12  636  560 (S25E32 L=250)   851    376        0 ( 568/ 421- 690/ 707)  45.0 E41-E24 S33-S15
     13  925  324 (S45E19 L=263)   466    109        0 ( 886/ 297- 945/ 367)  26.6 E23-E15 S47-S41
     14  621  342 (S39E42 L=240)   381     57        1 ( 591/ 326- 652/ 360) -26.6 E44-E39 S40-S37
     15  625  249 (S47E49 L=233)   567     97        1 ( 600/ 224- 648/ 283)  29.7 E51-E44 S48-S44
     16  741  257 (S48E38 L=244)   457     50        0 ( 715/ 235- 760/ 282)  21.8 E40-E35 S49-S45
  

    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