Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1075 1392 (N15E06 L=238)  1145    274        0 (1024/1270-1123/1529) -51.8 E12-E01 N08-N22
      2  371 1451 (N03E54 L=191)   319     75        0 ( 350/1410- 389/1491) -80.5 E55-E52 N00-N06
      3  460 1405 (N02E45 L=200)   936    273        0 ( 345/1381- 565/1425)  10.0 E53-E36 N00-N04
      4 1680 1327 (N29W34 L=279)   332    117        0 (1634/1311-1713/1339)  33.7 W30-W37 N26-N30
      5 1648 1247 (N22W33 L=277)   395     80        1 (1624/1196-1662/1286) -59.0 W30-W34 N19-N24
      6 1222 1175 (N06W07 L=252)   615    122        2 (1149/1124-1280/1220) -12.5 W02-W11 N04-N07
      7  432 1134 (S14E39 L=205)   438     76        0 ( 394/1107- 469/1166)  38.7 E41-E36 S15-S11
      8 1650  994 (N08W38 L=283)  1186    241        1 (1613/ 875-1685/1117) -45.0 W32-W44 N02-N13
      9  280 1078 (S20E50 L=194)   298    136        0 ( 247/1059- 321/1116)   8.1 E54-E47 S20-S17
     10  812  798 (S25E07 L=238)   850    282        0 ( 665/ 765- 945/ 846)  32.0 E16-W00 S29-S19
     11  349  754 (S38E41 L=203)  1482    388        2 ( 229/ 707- 466/ 785)  19.2 E53-E30 S43-S35
  

    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