Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1253 1610 (N36E04 L= 83)  1805    474        2 (1182/1406-1353/1810) -46.2 E16-W08 N25-N51
      2 1035 1541 (N24E14 L= 72)  1351    417        3 ( 901/1383-1174/1639) -16.3 E25-E01 N18-N28
      3  600 1565 (N15E42 L= 44)   277     99        0 ( 563/1528- 627/1587)  -9.5 E44-E38 N12-N15
      4  349 1472 (N04E58 L= 29)   362     80        0 ( 317/1445- 383/1489)  31.0 E59-E54 N01-N04
      5 1489 1368 (N27W18 L=105)   250     69        0 (1467/1335-1516/1396) -26.6 W16-W21 N26-N29
      6 1595 1287 (N26W28 L=114)   906    210        2 (1524/1241-1638/1370)  68.2 W22-W30 N22-N32
      7 1860 1306 (N38W57 L=143)  1108    122        1 (1826/1254-1905/1339)  26.6 W50-W68 N32-N41
      8  289 1285 (S08E55 L= 31)   357    111        1 ( 269/1234- 302/1330) -63.4 E58-E53 S11-S05
      9  717 1218 (S02E22 L= 64)   285     87        0 ( 682/1191- 751/1244)  26.6 E24-E20 S04-N00
     10 1518 1210 (N19W24 L=110)   336     71        1 (1472/1185-1553/1225)  26.6 W20-W26 N17-N20
     11  796 1170 (S03E16 L= 70)   232     33        1 ( 776/1145- 812/1197) -71.6 E17-E15 S04-S01
     12  762 1020 (S12E15 L= 72)   422    170        1 ( 681/ 982- 810/1070)  49.4 E19-E12 S15-S08
     13  406 1025 (S21E39 L= 48)   274     66        0 ( 377/1008- 444/1044)   0.0 E41-E35 S21-S19
     14 1434  893 (S01W26 L=113)  1179    275        3 (1382/ 773-1476/1026) -50.2 W21-W32 S06-N06
     15  380  896 (S29E39 L= 47)   280     70        1 ( 368/ 861- 394/ 934) -53.1 E40-E37 S30-S26
     16 1738  783 (N02W52 L=138)   905    280        1 (1703/ 637-1759/ 881) -42.5 W45-W59 S04-N07
  

    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