Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1154 1879 (N54E29 L= 42)   525     54        2 (1131/1855-1184/1899)   9.5 E33-E26 N51-N56
      2  644 1789 (N29E56 L= 15)   537     75        1 ( 618/1764- 679/1816) -10.3 E61-E50 N27-N30
      3  588 1663 (N19E49 L= 22)  1094    189        3 ( 554/1578- 653/1733) -24.0 E57-E39 N14-N23
      4 1236 1358 (N18W02 L= 73)   220     98        1 (1187/1347-1289/1370)  31.0 W00-W05 N15-N19
      5 1199 1307 (N14W02 L= 73)   246     27        0 (1175/1287-1222/1325) -18.4 W00-W03 N12-N14
      6 1231 1221 (N10W06 L= 77)   253     75        2 (1211/1174-1247/1253) -53.1 W04-W07 N07-N11
      7  178 1111 (S22E65 L=  6)   603     54        1 ( 157/1082- 198/1142)  63.4 E67-E62 S23-S19
      8 1621 1106 (N15W33 L=104)   242     45        2 (1605/1078-1638/1134) -45.0 W31-W34 N14-N16
      9  210 1044 (S25E59 L= 12)  1708    181        6 ( 164/ 964- 249/1126)  68.2 E65-E54 S30-S20
     10  376 1008 (S24E41 L= 30)   864    304        0 ( 259/ 925- 488/1056)  34.4 E52-E33 S31-S18
     11 1079  990 (S07W04 L= 75)   294    111        1 (1054/ 951-1121/1022)  63.4 W02-W06 S09-S05
  

    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