Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1771 (N49E05 L= 50)   508     87        0 (1248/1734-1313/1798) -29.7 E08-W00 N46-N51
      2  794 1713 (N29E35 L= 20)   321     49        0 ( 775/1680- 810/1743) -38.7 E37-E32 N26-N31
      3  638 1663 (N22E43 L= 12)   309     74        0 ( 623/1632- 659/1706) -39.8 E46-E40 N20-N25
      4 1602 1655 (N52W29 L= 84)   446     34        0 (1582/1637-1629/1673)  29.7 W26-W33 N49-N54
      5 1268 1586 (N33W01 L= 56)   406     92        0 (1228/1547-1305/1626) -26.6 E03-W03 N30-N35
      6 1324 1482 (N27W07 L= 62)   251     68        1 (1308/1446-1339/1519) -68.2 W05-W08 N24-N29
      7  194 1040 (S22E61 L=354)  1976    205        2 ( 141/ 962- 269/1097)  22.4 E69-E52 S27-S18
      8 1349  395 (S34W43 L= 99)   479     65        2 (1319/ 371-1390/ 419)   9.5 W39-W46 S35-S32
      9 1159  139 (S57W64 L=120)   825     84        0 (1144/  95-1173/ 179)  -9.5 W55-W85 S59-S53
  

    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