Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1772 (N33E42 L= 21)  2532    518        3 ( 750/1649- 829/1913) -39.8 E56-E31 N25-N45
      2 1498 1568 (N43W15 L= 78)   287     52        0 (1484/1538-1515/1591) -78.7 W13-W16 N39-N44
      3 1732 1399 (N39W39 L=101)  2669    378        5 (1627/1286-1801/1482) -25.0 W28-W46 N32-N43
      4  880 1242 (N03E14 L= 48)   196    101        0 ( 824/1230- 925/1258)  31.0 E17-E11 N01-N04
      5  446 1205 (S10E40 L= 22)   274     80        0 ( 409/1172- 478/1245)  63.4 E41-E38 S12-S06
      6 1429 1108 (N11W21 L= 83)   390    100        1 (1406/1051-1448/1161) -63.4 W18-W21 N07-N13
      7 1537 1117 (N14W27 L= 90)   651    211        1 (1485/1049-1595/1161)  76.0 W24-W31 N09-N17
      8  478  909 (S26E31 L= 31)   460    143        1 ( 419/ 891- 543/ 924)  18.4 E35-E26 S27-S23
      9  711  788 (S27E12 L= 50)   279     55        1 ( 683/ 767- 745/ 809)  45.0 E13-E10 S28-S24
     10  511  688 (S39E25 L= 38)   459     90        2 ( 474/ 669- 558/ 714)   7.1 E28-E20 S40-S37
     11 1354  285 (S36W51 L=114)   493     73        0 (1313/ 268-1399/ 302)  45.0 W48-W53 S38-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