Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1071 1776 (N42E25 L=220)  1449    273        1 (1018/1655-1151/1847) -69.8 E29-E18 N32-N51
      2 1416 1482 (N32W10 L=256)   797    128        1 (1386/1411-1434/1542) -78.7 W08-W12 N27-N37
      3 1634 1403 (N35W29 L=275)  1381    310        0 (1560/1260-1693/1497) -35.5 W21-W36 N28-N41
      4 1068  581 (S29W17 L=262)   207     54        0 (1054/ 552-1082/ 610) -71.6 W15-W18 S30-S27
      5  847  548 (S38W04 L=250)   308     90        0 ( 812/ 511- 881/ 582)  74.1 W02-W06 S41-S34
      6 1339  473 (S26W39 L=284)   258     98        0 (1299/ 459-1385/ 489)  38.7 W35-W40 S27-S23
      7 1015  419 (S41W22 L=267)   471    101        0 ( 971/ 388-1046/ 461)  81.9 W19-W23 S44-S37
      8 1234  286 (S41W47 L=292)   485    100        0 (1179/ 258-1272/ 310)  11.3 W41-W52 S42-S39
      9 1476  275 (S32W68 L=313)  1190    116        2 (1411/ 249-1527/ 300)  12.5 W60-W79 S34-S29
  

    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