Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1097 (S16E44 L=207)  1224    290        3 ( 281/1007- 452/1205)  59.0 E49-E39 S22-S07
      2  260 1131 (S16E54 L=197)   323     69        1 ( 245/1097- 282/1161)  76.0 E55-E52 S17-S13
      3 1873 1052 (N21W58 L=310)   460    123        0 (1835/1029-1893/1085)  45.0 W54-W60 N19-N24
      4  182  812 (S37E58 L=193)  1059    268        1 ( 134/ 775- 222/ 858)  31.0 E63-E53 S39-S33
      5  411  778 (S34E35 L=216)   697    149        0 ( 354/ 739- 473/ 822) -10.3 E40-E29 S35-S32
      6  811  759 (S25E06 L=246)  1657    446        3 ( 574/ 677-1002/ 809)   9.8 E22-W07 S28-S21
      7 1412  659 (S14W33 L=285)  1537    371        3 (1247/ 627-1586/ 696)  18.4 W21-W45 S17-S07
      8  572  537 (S46E19 L=233)   729    120        2 ( 531/ 507- 628/ 583)  -6.3 E22-E13 S48-S42
      9 1251  381 (S35W37 L=288)   304     27        1 (1229/ 362-1276/ 398)  71.6 W34-W38 S36-S33
     10 1403  277 (S36W54 L=306)   438     96        0 (1355/ 266-1440/ 294)  38.7 W51-W57 S38-S34
  

    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