Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  577 1664 (N19E49 L=225)   536     79        1 ( 541/1637- 618/1695)  45.0 E51-E47 N17-N21
      2 1224 1632 (N35E07 L=267)   288     39        1 (1207/1604-1243/1654) -45.0 E08-E05 N33-N37
      3  358 1412 (S01E54 L=219)   530     83        0 ( 344/1365- 375/1450) -76.0 E56-E52 S03-N01
      4  829 1373 (N08E21 L=253)   570    184        1 ( 759/1319- 905/1442) -14.0 E26-E14 N06-N10
      5 1298 1217 (N12W10 L=283)  1294    322        0 (1261/1036-1322/1362) -57.1 W04-W15 N01-N19
      6  214 1196 (S16E60 L=214)   770    186        3 ( 182/1159- 273/1235)  45.0 E63-E55 S18-S12
      7  127 1196 (S17E71 L=202)   887     95        1 ( 109/1163- 144/1229)  90.0 E73-E68 S18-S14
      8 1504  940 (N03W29 L=303)   376    116        0 (1466/ 891-1524/1006) -80.5 W27-W30 N00-N06
      9  372  935 (S28E40 L=234)  1263    368        1 ( 245/ 888- 502/1001)  23.2 E51-E30 S32-S21
     10  995  882 (S15W02 L=276)  1096    349        1 ( 810/ 837-1154/ 923)  16.7 E09-W11 S18-S11
     11  208  682 (S48E56 L=217)   532     91        0 ( 187/ 666- 239/ 705)   9.5 E58-E52 S49-S46
     12 1433  345 (S30W52 L=326)   350    101        0 (1377/ 323-1466/ 371)  71.6 W49-W54 S33-S27
  

    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