Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  539 1721 (N47E51 L=129)  3654    512        6 ( 422/1643- 668/1773)  16.7 E60-E36 N40-N51
      2 1111 1645 (N45W05 L=184)   282     59        1 (1094/1611-1130/1669)  63.4 W03-W06 N42-N46
      3  822 1632 (N42E20 L=160)   394     72        0 ( 781/1610- 851/1661)  45.0 E22-E17 N40-N44
      4 1225 1550 (N38W13 L=193)   305     66        1 (1197/1532-1259/1579) -21.8 W11-W16 N36-N39
      5  418 1362 (N21E45 L=135)   428     89        1 ( 391/1336- 446/1392) -26.6 E47-E41 N19-N23
      6  356 1270 (N15E48 L=131)   988    219        3 ( 288/1235- 430/1295)   4.4 E55-E42 N12-N16
      7  320 1189 (N10E51 L=129)  1167    298        1 ( 246/1141- 422/1244) -16.4 E58-E41 N06-N12
      8  926 1060 (N05E06 L=173)   235     74        0 ( 900/1028- 964/1095) -36.9 E07-E03 N02-N06
      9  600  939 (S04E27 L=153)   240     82        0 ( 582/ 899- 621/ 985)  68.2 E27-E25 S06-S01
     10  370  818 (S13E45 L=135)  1468    224        5 ( 324/ 736- 410/ 915) -85.2 E48-E41 S18-S06
     11 1854  843 (S07W65 L=244)   711    102        1 (1829/ 804-1871/ 880) -51.3 W60-W67 S09-S04
     12  948  780 (S13E04 L=176)   283     57        1 ( 922/ 743- 976/ 816)  63.4 E05-E02 S14-S10
     13 1409  583 (S24W28 L=208)   303     77        0 (1373/ 569-1451/ 599)  -9.5 W25-W31 S25-S23
     14 1218  432 (S36W17 L=197)   856    163        2 (1144/ 410-1295/ 464)  -9.5 W11-W23 S37-S33
     15  808  403 (S40E15 L=165)  1371    321        5 ( 673/ 357- 943/ 442)  12.3 E26-E03 S44-S36
     16 1425  380 (S40W36 L=216)   346     46        0 (1404/ 353-1446/ 402)  76.0 W34-W37 S41-S37
  

    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