Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1154 1696 (N47E02 L=118)  1138    247        1 (1121/1621-1205/1775) -85.9 E05-W01 N41-N55
      2  569 1579 (N28E42 L= 78)  3004    513        7 ( 421/1468- 665/1750)  79.2 E51-E37 N19-N41
      3  389 1660 (N31E64 L= 56)   873    171        0 ( 376/1629- 405/1699) -59.0 E66-E60 N28-N33
      4 1509 1579 (N44W30 L=151)   876    171        2 (1447/1550-1564/1611)  22.6 W24-W36 N41-N47
      5  578 1403 (N17E35 L= 85)   487    170        0 ( 517/1368- 621/1459)  50.2 E39-E32 N14-N20
      6 1502 1419 (N32W28 L=148)   477     98        2 (1479/1376-1544/1456)  80.5 W26-W31 N28-N34
      7  387 1361 (N11E49 L= 71)   436    116        3 ( 366/1320- 411/1402) -33.7 E52-E45 N09-N13
      8  724 1248 (N09E22 L= 98)   306     80        1 ( 702/1212- 751/1285) -45.0 E23-E19 N07-N11
      9 1216 1212 (N14W09 L=129)   534    121        1 (1170/1173-1257/1268)  90.0 W06-W11 N11-N17
     10 1242 1083 (N07W12 L=132)   854    343        1 (1108/ 955-1334/1174) -29.4 W03-W19 N00-N10
     11  252  950 (S15E54 L= 66)   818    223        0 ( 215/ 865- 270/1058)  77.0 E57-E51 S21-S08
     12 1665  913 (N02W44 L=164)   551    104        1 (1646/ 865-1684/ 973) -63.4 W42-W46 N00-N06
     13  397  917 (S15E40 L= 80)   312     87        1 ( 378/ 886- 438/ 946)  36.9 E42-E37 S17-S13
     14  841  883 (S11E09 L=111)   334    158        0 ( 758/ 864- 916/ 905)  11.3 E14-E04 S12-S08
     15 1218  770 (S12W16 L=136)   784     94        4 (1163/ 710-1260/ 829) -39.8 W11-W18 S15-S08
     16 1742  684 (S10W57 L=177)   488     82        0 (1720/ 647-1762/ 728)  78.7 W54-W58 S12-S07
     17  701  555 (S34E16 L=104)   263     34        1 ( 681/ 536- 720/ 573) -33.7 E17-E14 S35-S32
  

    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