Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1172 1734 (N42W36 L=132)   608    158        0 (1103/1702-1240/1764)  -9.5 W30-W42 N40-N42
      2 1556 1486 (N15W47 L=143)  1445    362        2 (1415/1433-1713/1540) -39.8 W39-W57 N07-N22
      3 1155 1482 (N27W21 L=117)   328    127        0 (1095/1460-1213/1507)  -7.1 W17-W25 N25-N28
      4 1398 1389 (N15W32 L=128)   321     57        1 (1365/1367-1430/1417) -53.1 W30-W33 N12-N16
      5 1880 1141 (S13W62 L=158)  2060    406        7 (1837/ 997-1917/1256)  54.5 W54-W69 S22-S07
      6 1243 1066 (N02W13 L=109)   268     97        1 (1183/1047-1286/1094) -38.7 W10-W15 N00-N04
      7  365  969 (N18E43 L= 53)   763    203        3 ( 310/ 924- 456/1017) -35.0 E46-E36 N14-N21
      8 1690  958 (S17W39 L=135)   676    141        1 (1641/ 906-1737/1006) -58.0 W35-W41 S20-S12
      9 1720  830 (S26W39 L=135)   314     52        1 (1703/ 805-1738/ 853)  33.7 W37-W40 S27-S24
     10 1374  773 (S18W13 L=109)   241    100        0 (1319/ 746-1421/ 805)   0.0 W09-W16 S18-S17
     11 1203  671 (S19W01 L= 97)   349    145        0 (1174/ 634-1238/ 709)  26.6 E01-W03 S21-S16
     12 1589  652 (S33W26 L=122)   276     46        1 (1564/ 628-1610/ 676)  26.6 W23-W27 S33-S30
     13 1353  513 (S34W06 L=102)   987    172        0 (1310/ 446-1389/ 590)  45.0 W01-W09 S37-S29
     14  795  415 (S24E32 L= 64)   933    224        2 ( 765/ 328- 841/ 501)  57.5 E35-E26 S31-S19
     15 1531  412 (S49W18 L=114)  2028    259        3 (1457/ 351-1605/ 487)  38.7 W11-W26 S54-S43
     16 1442  366 (S49W09 L=105)   855    124        3 (1399/ 328-1483/ 406) -32.0 W05-W13 S52-S44
     17  528  329 (S23E56 L= 40)  1759    282        2 ( 501/ 261- 569/ 405)  47.7 E63-E52 S28-S17
  

    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