Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  498 1338 (N33E27 L=178)   999    230        5 ( 414/1289- 564/1375) -26.6 E34-E22 N29-N36
      2  142  926 (N15E66 L=139)   510    104        1 ( 124/ 894- 170/ 956)  26.6 E69-E61 N13-N16
      3  670  682 (S08E29 L=177)  3499    765       10 ( 496/ 458- 811/ 845) -17.1 E40-E19 S21-N01
      4 1104  697 (S17E03 L=202)   189     49        3 (1074/ 678-1129/ 722) -18.4 E04-E01 S18-S15
      5  527  591 (S11E41 L=164)   287    127        2 ( 501/ 544- 562/ 627)  29.7 E44-E37 S13-S09
      6  886  498 (S24E22 L=184)   518    148        2 ( 830/ 461- 969/ 535)   0.0 E25-E15 S25-S22
      7  493  453 (S18E50 L=155)   935    232        2 ( 432/ 383- 551/ 524) -78.7 E53-E46 S23-S13
  

    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