Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  699 1628 (N19E38 L= 81)   683    171        0 ( 674/1545- 736/1694) -63.4 E40-E34 N14-N24
      2 1064 1615 (N29E16 L=104)   467    126        2 (1029/1565-1091/1673) -36.9 E20-E12 N25-N31
      3  848 1531 (N17E25 L= 94)   884    272        0 ( 766/1423- 899/1622) -45.0 E32-E20 N10-N21
      4  323 1436 (N00E59 L= 61)   938    266        0 ( 263/1384- 373/1469)  -3.6 E67-E51 S02-N01
      5  603 1322 (S01E32 L= 87)   266     65        0 ( 577/1294- 632/1353) -11.3 E34-E29 S01-N00
      6  648 1214 (S06E26 L= 93)   384     79        1 ( 617/1166- 679/1257)  80.5 E27-E24 S09-S03
      7  206 1033 (S26E57 L= 63)  3097    437        7 ( 120/ 880- 266/1196)  65.8 E64-E52 S35-S15
      8 1385 1101 (N08W18 L=137)   568    178        0 (1306/1078-1470/1127)  15.3 W12-W23 N05-N09
      9 1383 1056 (N05W19 L=138)   210     58        1 (1351/1029-1409/1086) -14.0 W16-W20 N03-N06
     10  432 1043 (S20E37 L= 82)   238     74        1 ( 416/1002- 445/1083)  80.5 E38-E35 S23-S17
     11  653 1019 (S17E21 L= 98)   403    107        1 ( 605/ 975- 701/1051)  -8.1 E24-E17 S18-S15
     12 1595  970 (N07W35 L=154)   488    131        0 (1535/ 929-1643/1009)  49.4 W31-W37 N03-N10
     13  740  905 (S21E13 L=106)   623    164        1 ( 670/ 846- 813/ 961)  59.0 E16-E09 S25-S15
     14 1824  772 (N05W61 L=181)   457    129        0 (1788/ 752-1857/ 795)   6.3 W56-W65 N03-N06
     15 1108  555 (S30W20 L=139)   798    177        1 (1023/ 522-1195/ 602)  45.0 W15-W24 S34-S25
     16 1713  424 (S16W72 L=191)  2283    281        3 (1666/ 393-1770/ 458)  38.7 W66-W78 S19-S11
     17 1213  271 (S43W44 L=163)   577     75        1 (1181/ 248-1242/ 305)  21.8 W41-W47 S44-S40
  

    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