Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  893 1701 (N47E17 L=252)   350     64        0 ( 867/1683- 933/1716)  -9.5 E19-E13 N45-N48
      2 1012 1657 (N44E06 L=264)   496    124        0 ( 954/1638-1064/1679)  24.0 E10-E01 N42-N46
      3 1281 1612 (N42W17 L=287)  1979    315        2 (1121/1576-1426/1643)  -4.4 W03-W29 N39-N44
      4 1732 1571 (N41W68 L=338)   877     61        0 (1696/1557-1770/1585)   0.0 W61-W81 N39-N41
      5  342 1525 (N29E61 L=208)   438     46        1 ( 327/1501- 360/1551) -26.6 E64-E58 N27-N30
      6 1699 1507 (N36W56 L=326)   995     98        1 (1670/1452-1740/1545)  26.6 W51-W65 N32-N39
      7 1813 1325 (N24W62 L=332)   506    115        0 (1788/1268-1830/1382) -81.9 W61-W63 N20-N27
      8 1764 1303 (N23W55 L=325)   698    122        0 (1732/1249-1791/1352) -63.4 W52-W57 N19-N25
      9 1139 1160 (N11W06 L=276)  1806    423        2 (1087/ 946-1166/1341) -85.0 W02-W08 S01-N22
     10  433 1148 (N06E40 L=230)  1686    344        3 ( 383/1001- 473/1308) -70.6 E46-E36 S02-N15
     11  693 1258 (N15E23 L=247)   247     85        0 ( 650/1240- 733/1277)  21.8 E25-E20 N13-N16
     12  604 1183 (N10E28 L=242)   638    176        0 ( 577/1112- 635/1254) -83.7 E29-E25 N05-N14
     13  672  763 (S16E21 L=248)  1204    403        0 ( 607/ 607- 712/ 922)  74.1 E26-E18 S26-S05
     14 1646  451 (S32W57 L=327)  1695    103        3 (1590/ 417-1705/ 481)  -6.3 W48-W66 S34-S29
     15 1455  369 (S39W42 L=311)   557     65        0 (1418/ 351-1493/ 386)   0.0 W37-W45 S40-S37
     16  675  285 (S53E30 L=240)   484     42        1 ( 656/ 261- 698/ 312)  63.4 E33-E27 S55-S49
     17  588  278 (S55E43 L=227)  1031     96        1 ( 536/ 253- 633/ 308)  15.3 E57-E34 S58-S50
  

    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