Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  867 1702 (N41E26 L=124)   277     41        2 ( 845/1682- 888/1721)  14.0 E28-E24 N39-N42
      2 1258 1596 (N41W07 L=158)   810    186        1 (1189/1530-1316/1639)  45.0 W02-W11 N34-N44
      3  545 1404 (N16E38 L=113)   378    144        4 ( 496/1378- 589/1443)  59.0 E40-E34 N13-N18
      4  732 1284 (N11E22 L=128)   897    219        3 ( 675/1217- 795/1339)  66.0 E25-E18 N06-N15
      5  139 1014 (S14E69 L= 82)   840    140        3 ( 115/ 967- 158/1049) -78.7 E72-E66 S16-S11
      6  435 1014 (S09E38 L=112)   285     98        0 ( 391/ 992- 472/1043)  26.6 E41-E35 S10-S07
      7  747  924 (S10E15 L=135)   929    213        2 ( 712/ 826- 790/1018) -70.0 E18-E11 S15-S04
      8 1236  992 (N01W13 L=164)   250     46        0 (1205/ 974-1264/1010) -14.0 W11-W15 N00-N01
      9 1434  912 (S01W27 L=178)   434     96        1 (1391/ 885-1485/ 935)  26.6 W24-W30 S02-N01
     10  294  595 (S39E52 L= 98)   687    177        1 ( 261/ 572- 333/ 623)  29.1 E56-E47 S41-S36
     11 1311  427 (S33W32 L=182)   325     66        0 (1277/ 406-1351/ 446)  36.9 W30-W34 S35-S31
     12 1414  377 (S35W43 L=193)   620    134        1 (1373/ 355-1468/ 410)   0.0 W38-W47 S36-S32
  

    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