Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  400 1411 (N15E50 L=331)   862    202        3 ( 380/1322- 437/1508) -53.1 E56-E45 N09-N21
      2 1644 1462 (N37W43 L= 65)   348     54        1 (1620/1446-1671/1485) -21.8 W40-W46 N36-N38
      3 1543 1371 (N29W32 L= 53)   400    121        0 (1497/1333-1583/1411) -29.7 W28-W35 N27-N31
      4  599 1259 (N09E30 L=351)   716    221        1 ( 547/1142- 638/1340)  81.3 E31-E28 N01-N14
      5 1826 1233 (N24W60 L= 81)  1113    186        2 (1754/1163-1882/1290) -20.6 W51-W67 N20-N26
      6  312 1155 (S01E51 L=331)   533     86        0 ( 285/1130- 347/1188) -18.4 E53-E47 S02-N00
      7  811  697 (S23E09 L= 12)  1877    628        2 ( 578/ 598-1070/ 804)  33.7 E25-W05 S32-S12
      8  352  665 (S32E47 L=335)   415    101        1 ( 303/ 654- 402/ 676)   5.2 E52-E41 S32-S30
      9 1226  581 (S24W20 L= 41)   393    102        0 (1175/ 567-1275/ 596)  23.2 W16-W23 S25-S22
     10 1083  453 (S35W14 L= 35)   295     91        0 (1039/ 436-1128/ 469)   0.0 W10-W17 S36-S34
  

    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