Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  872 1761 (N45E31 L=209)   596     50        1 ( 849/1732- 903/1792) -45.0 E33-E27 N42-N47
      2  283 1519 (N17E70 L=170)  1042     78        1 ( 254/1488- 321/1553)  -6.0 E81-E62 N15-N18
      3 1752 1487 (N42W56 L=296)   597     84        0 (1713/1460-1799/1506)  -4.4 W50-W64 N41-N43
      4 1274 1323 (N21W11 L=250)   280     94        1 (1237/1281-1306/1354) -31.0 W07-W13 N19-N23
      5 1078 1208 (N11W00 L=240)  1547    471        4 ( 997/1117-1186/1278)  45.0 E04-W05 N05-N16
      6  609 1223 (N04E29 L=211)   321    185        1 ( 541/1193- 688/1251)   0.0 E33-E23 N03-N06
      7  829 1180 (N05E14 L=226)   376    135        0 ( 816/1112- 843/1245) -69.4 E15-E12 N01-N09
      8 1342 1070 (N07W18 L=258)   297    120        0 (1295/1040-1395/1089)   0.0 W15-W22 N05-N08
      9  723  995 (S07E18 L=222)   218     56        0 ( 702/ 964- 739/1016) -33.7 E19-E16 S08-S05
     10 1079  829 (S11W06 L=246)  1087    373        0 (1009/ 686-1172/ 981)  81.5 W03-W09 S20-N00
     11  188  895 (S22E62 L=178)  1238    174        0 ( 124/ 842- 252/ 943)  22.2 E75-E53 S26-S17
     12 1510  652 (S15W38 L=278)   380     90        1 (1465/ 620-1543/ 680) -15.9 W34-W41 S16-S13
     13  523  368 (S55E35 L=205)  1925    245        5 ( 438/ 303- 610/ 455) -32.0 E42-E25 S59-S48
  

    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