Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1464 1731 (N32W56 L=246)   962    191        1 (1401/1707-1529/1751) -26.6 W50-W61 N29-N34
      2 1525 1437 (N14W43 L=233)  2801    665        5 (1263/1399-1743/1486) -23.2 W25-W61 N05-N21
      3 1071 1455 (N29W16 L=206)   280     77        0 (1035/1424-1115/1477)   0.0 W12-W19 N28-N30
      4  763 1281 (N28E09 L=181)   491    134        1 ( 728/1214- 805/1342)  83.7 E10-E05 N23-N32
      5 1882 1208 (S10W65 L=255)   587    141        1 (1858/1148-1897/1246)  45.0 W62-W68 S13-S07
      6 1185 1174 (N11W13 L=203)   342    138        0 (1171/1100-1212/1235)  49.4 W10-W16 N06-N13
      7 1809 1106 (S12W53 L=243)  1397    326        6 (1756/ 989-1868/1169) -69.0 W48-W59 S20-S07
      8 1223 1045 (N03W11 L=201)   354     90        1 (1177/1025-1273/1071) -45.0 W09-W13 N00-N04
      9  359  960 (N19E44 L=146)   361     78        1 ( 349/ 914- 372/1005)  78.7 E45-E42 N16-N21
     10 1582  803 (S22W27 L=217)  1780    349        2 (1513/ 659-1619/ 960)  84.3 W23-W29 S31-S10
     11 1392  676 (S23W12 L=202)   385    109        0 (1367/ 616-1430/ 727)  39.8 W09-W15 S26-S21
     12  530  544 (S09E44 L=146)   345    128        0 ( 477/ 523- 580/ 586) -39.8 E46-E40 S10-S05
     13  297  537 (S05E67 L=123)  1204    156        1 ( 264/ 503- 349/ 581) -39.8 E71-E61 S07-S02
     14  728  469 (S18E34 L=156)   498    195        1 ( 686/ 401- 762/ 539)  28.6 E39-E28 S21-S14
     15  887  324 (S32E31 L=159)   368     70        0 ( 851/ 305- 913/ 349)   9.5 E34-E28 S32-S30
     16  737  241 (S34E47 L=143)   798    128        2 ( 716/ 203- 763/ 283)  33.7 E51-E42 S36-S30
  

    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