Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  900 1526 (N30E17 L=243)   299     87        1 ( 875/1483- 927/1568) -39.8 E19-E13 N27-N32
      2  656 1458 (N21E31 L=228)   697    119        2 ( 631/1400- 693/1509)  90.0 E33-E29 N17-N25
      3  432 1472 (N19E49 L=211)   436    115        1 ( 402/1428- 451/1507) -39.8 E52-E46 N16-N21
      4  965 1381 (N21E09 L=251)   280     90        1 ( 933/1336- 981/1414) -45.0 E11-E07 N18-N23
      5  548 1312 (N11E35 L=225)   733    153        0 ( 513/1253- 591/1384) -41.2 E39-E31 N07-N14
      6 1042 1223 (N12E02 L=258)   546    182        3 ( 941/1200-1127/1255)   0.0 E08-W03 N11-N13
      7 1147  981 (S01W08 L=268)  1233    317        1 (1080/ 853-1236/1117) -63.4 W03-W15 S07-N08
      8  263 1024 (S11E54 L=206)   777    163        3 ( 247/ 938- 280/1105) -78.7 E55-E51 S15-S05
      9  604  926 (S12E25 L=235)   505    114        2 ( 554/ 902- 673/ 942)   7.1 E28-E20 S12-S09
     10  891  875 (S11E06 L=254)   467    180        1 ( 847/ 795- 970/ 923)  49.4 E08-E01 S15-S07
     11 1223  707 (S16W17 L=277)   909    208        1 (1150/ 601-1288/ 787)  69.0 W14-W19 S23-S10
     12 1049  538 (S30W09 L=269)   493    161        0 ( 991/ 515-1122/ 566)  24.0 W05-W14 S32-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