Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1067 1772 (N45E07 L=258)  1452    256        1 ( 930/1732-1177/1811)  -5.2 E20-W02 N43-N47
      2 1547 1629 (N41W35 L=300)   476     79        0 (1506/1604-1588/1650)  -7.1 W30-W38 N39-N42
      3  804 1475 (N20E19 L=247)   930    240        1 ( 727/1389- 883/1547) -31.6 E25-E12 N15-N24
      4  244 1495 (N18E69 L=196)  1035     78        0 ( 225/1466- 276/1537) -24.4 E76-E62 N16-N21
      5 1496 1508 (N31W28 L=293)   254     68        0 (1477/1479-1524/1527)  45.0 W25-W30 N28-N32
      6 1347 1503 (N28W16 L=282)   334     70        1 (1308/1480-1379/1525)   0.0 W13-W18 N26-N29
      7 1273 1463 (N24W11 L=277)   274     82        1 (1240/1437-1312/1493) -21.8 W08-W13 N22-N26
      8 1502 1430 (N25W27 L=293)   375     43        0 (1478/1399-1530/1460) -53.1 W25-W29 N23-N27
      9 1684 1369 (N24W42 L=308)   314     56        1 (1660/1342-1713/1388)  31.0 W40-W45 N22-N25
     10  960 1348 (N13E07 L=258)   190     29        0 ( 933/1329- 982/1363)   0.0 E09-E05 N12-N14
     11 1548 1282 (N16W31 L=296)   413    118        0 (1529/1225-1578/1339)  69.4 W29-W33 N12-N20
     12 1697 1225 (N15W43 L=308)   685    178        0 (1651/1148-1742/1300) -48.8 W38-W46 N11-N19
     13 1885  801 (S06W69 L=334)   624     85        0 (1868/ 764-1904/ 844) -45.0 W65-W73 S08-S03
     14  233  742 (S28E60 L=205)  1104    148        0 ( 209/ 675- 274/ 808)  33.7 E66-E53 S32-S24
     15 1572  746 (S15W39 L=305)   344    108        0 (1517/ 733-1629/ 761)   5.7 W34-W44 S15-S14
     16 1761  673 (S16W59 L=324)  1924    231        6 (1718/ 593-1832/ 735) -66.0 W54-W68 S21-S12
     17  991  632 (S30W02 L=268)  1300    350        3 ( 870/ 519-1118/ 734)  50.9 E05-W09 S38-S21
  

    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