Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1416 1459 (N18W36 L=326)   263    101        2 (1370/1436-1449/1490)   0.0 W33-W39 N16-N19
      2 1119 1374 (N22W16 L=306)   194     61        0 (1096/1348-1150/1393)   0.0 W13-W17 N21-N22
      3 1766 1377 (N03W60 L=350)   574    239        2 (1703/1361-1827/1393) -15.9 W54-W68 N00-N05
      4 1228 1141 (N06W14 L=304)   782    249        2 (1104/1070-1358/1171) -25.0 W07-W22 N00-N10
      5  819 1071 (N13E10 L=280)   403    169        1 ( 736/1020- 893/1109)   5.2 E16-E05 N12-N14
      6  608  981 (N13E25 L=265)   202    102        1 ( 550/ 968- 650/ 993) -33.7 E28-E22 N11-N15
      7  187  877 (N16E62 L=228)  2624    557        4 (  85/ 823- 301/ 933) -14.0 E75-E51 N11-N20
      8 1343  650 (S25W08 L=298)  1441    362        6 (1281/ 496-1404/ 792)  45.0 W00-W15 S33-S17
      9  695  545 (S13E33 L=257)   613    178        0 ( 639/ 493- 765/ 613) -71.6 E35-E30 S17-S08
     10  605  408 (S20E45 L=245)   322     70        2 ( 582/ 385- 633/ 429) -14.0 E47-E43 S21-S18
     11  587  352 (S23E50 L=240)   355     71        1 ( 560/ 336- 611/ 376) -56.3 E52-E47 S24-S21
  

    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