Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  331 1457 (N21E60 L=241)  1536    211        4 ( 294/1395- 380/1537) -66.0 E65-E54 N16-N25
      2  522 1302 (N13E36 L=265)   293     70        2 ( 495/1285- 557/1327)  11.3 E38-E33 N11-N14
      3  638 1060 (S00E25 L=277)   281     46        1 ( 610/1035- 667/1084)   0.0 E26-E22 S01-N00
      4  672  898 (S10E21 L=280)   311     99        0 ( 641/ 867- 718/ 916)  21.8 E22-E17 S11-S08
      5  532  895 (S12E31 L=271)   237     46        1 ( 506/ 873- 556/ 915)   0.0 E32-E28 S12-S10
      6  409  823 (S18E40 L=261)   542    143        1 ( 341/ 791- 455/ 862)  33.7 E46-E36 S20-S14
      7  692  783 (S17E19 L=282)   473    154        1 ( 624/ 763- 745/ 820) -14.0 E23-E15 S18-S15
      8 1380  719 (S14W26 L=328)   604    175        0 (1357/ 657-1412/ 790)  90.0 W24-W29 S17-S09
      9 1113  714 (S17W09 L=310)   383    176        0 (1058/ 686-1178/ 761)  35.5 W05-W12 S18-S13
     10  566  597 (S31E29 L=273)   777    246        2 ( 481/ 560- 640/ 655) -24.4 E35-E22 S33-S26
     11 1368  374 (S38W37 L=339)   891     81        3 (1332/ 346-1407/ 407)  59.0 W33-W40 S40-S35
  

    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