Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1691 (N51W17 L=290)   542    122        2 ( 887/1644- 969/1736)  24.4 W10-W21 N48-N54
      2 1271 1652 (N37W38 L=311)   273     30        0 (1258/1624-1287/1674)  45.0 W36-W41 N35-N38
      3 1360 1654 (N34W45 L=318)   310     51        0 (1333/1637-1393/1666) -26.6 W42-W48 N32-N35
      4 1648 1563 (N20W62 L=335)  2018    142        4 (1596/1515-1726/1608) -35.0 W57-W71 N16-N23
      5 1485 1555 (N25W47 L=320)   405     90        1 (1462/1508-1501/1598)  63.4 W44-W49 N21-N27
      6  399 1325 (N39E40 L=233)  1215    172        0 ( 319/1274- 472/1358) -25.2 E49-E32 N34-N42
      7 1021 1247 (N20W05 L=278)   280     92        1 ( 984/1216-1063/1272) -45.0 W03-W07 N17-N22
      8  819 1250 (N25E07 L=266)   781    220        3 ( 708/1218- 906/1270) -26.6 E14-E01 N21-N27
      9 1641 1101 (S04W40 L=313)   996    147        3 (1601/1027-1669/1174)  49.4 W35-W43 S08-N00
     10 1195 1128 (N09W12 L=285)   413    139        1 (1137/1095-1235/1167)   8.1 W08-W15 N07-N10
     11  691 1106 (N19E18 L=255)   353    114        2 ( 648/1052- 722/1155)  63.4 E21-E15 N16-N22
     12  568 1080 (N20E27 L=246)   239    140        0 ( 539/1052- 598/1114)  26.6 E29-E24 N18-N21
     13  386 1049 (N22E42 L=231)  1041    173        1 ( 322/ 988- 454/1104) -45.0 E47-E36 N17-N26
     14 1534  565 (S35W23 L=295)  5184    863        7 (1324/ 325-1737/ 901)  31.8 W01-W43 S49-S18
     15 1082  828 (S05E01 L=272)   225     65        0 (1063/ 801-1106/ 860) -76.0 E02-W00 S07-S03
     16 1078  563 (S21E08 L=265)   502    189        0 (1020/ 504-1106/ 634)  41.2 E13-E05 S23-S16
     17  928  324 (S34E27 L=246)   720    155        0 ( 873/ 280- 983/ 376)  21.0 E33-E20 S36-S31
     18 1159  231 (S50E17 L=256)   778    101        2 (1110/ 207-1210/ 254)  -6.3 E21-E11 S52-S47
  

    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