Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  855 1627 (N33E26 L=196)   757    143        0 ( 800/1583- 901/1655)   7.1 E29-E21 N30-N35
      2 1188 1519 (N33W00 L=223)  1026    188        2 (1118/1451-1257/1572) -24.4 E05-W06 N29-N35
      3 1344 1285 (N20W15 L=237)   321    102        0 (1326/1241-1364/1327) -68.2 W13-W16 N18-N23
      4 1114 1225 (N12W02 L=224)   557    219        1 (1023/1189-1217/1261)   0.0 E04-W08 N11-N13
      5 1603 1202 (N21W35 L=257)   326     68        1 (1571/1175-1636/1222) -21.8 W32-W37 N19-N22
      6 1746 1055 (N15W49 L=271)  1011    184        3 (1717/ 983-1798/1128) -48.8 W45-W54 N10-N19
      7  227  950 (S20E58 L=165)  3154    367        3 ( 117/ 824- 326/1084)  30.5 E83-E47 S30-S10
      8 1792  969 (N10W55 L=277)   816    212        3 (1736/ 881-1834/1062)  63.4 W50-W58 N04-N16
      9 1060  736 (S18W08 L=230)   335     42        1 (1033/ 710-1087/ 758) -26.6 W05-W09 S19-S16
     10 1590  676 (S11W45 L=267)  1012    179        0 (1504/ 647-1652/ 712)  20.0 W39-W50 S14-S08
     11  905  608 (S29W01 L=223)   336    137        0 ( 845/ 578- 967/ 632)  33.7 E02-W04 S31-S26
     12 1624  585 (S15W52 L=274)   362     87        0 (1605/ 555-1641/ 624) -63.4 W48-W53 S17-S13
     13  765  302 (S57E03 L=219)  1209    261        4 ( 622/ 281- 872/ 336)  15.9 E21-W08 S62-S52
  

    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