Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  325 1256 (N32E43 L=179)  1084    185        4 ( 267/1198- 378/1294) -37.9 E47-E37 N28-N35
      2  394  705 (S01E47 L=174)  2330    418        5 ( 313/ 597- 518/ 773)  -6.3 E56-E38 S06-N02
      3  848  687 (S12E18 L=203)   263    109        0 ( 793/ 674- 909/ 698) -15.9 E21-E14 S13-S10
      4  549  526 (S15E42 L=179)  2109    424        5 ( 468/ 398- 613/ 665)  37.9 E54-E35 S21-S06
      5  352  489 (S13E61 L=160)  1456    261        5 ( 299/ 426- 393/ 557)  41.6 E67-E53 S17-S09
      6  693  406 (S25E38 L=183)  1233    399        2 ( 617/ 289- 790/ 497)  25.5 E49-E28 S32-S21
      7  417  383 (S21E61 L=160)   520    219        1 ( 373/ 360- 458/ 418)   5.2 E66-E55 S22-S19
  

    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